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Abstract 

We investigate the effect of weak disorder on different coupled g-state Potts models 
with q < 4 using two loops renormalisation group. This study presents new examples 
of first order transitions driven by randomness. We found that weak disorder makes 
the models decouple. Therefore, it appears that no relations emerge, at a perturbation 
level, between the disordered q\ x g 2 -state Potts model and the two disordered q±, 
g2-state Potts models (qi ^ (72), despite their central charges are similar according 
to recent numerical investigations. Nevertheless, when two g-state Potts models are 
considered (q > 2), the system remains always driven in a strong coupling regime, 
violating apparently the Imry-Wortis argument. 
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The effect of weak randomness on continuous phase transitions has been much studied either 
analytically |l], 0, or numerically ||. For most of these cases, the Harris criterion j| pro- 
vides us a good method to see if disorder will change the universality class of the transition. 
However, it has been realized recently, that a weak bond randomness can also have strong 
effects on first order transitions 2D systems and can induce a second order phase transition 
in a system that would undergo a first order one ||. This result has been established more 
rigorously in ||. This has been first tested with success on systems presenting a fluctuation 
driven first order transition, namely transitions which are expected to be continuous at a 
mean field level, but become first order when transitions are incorporated. Hence, Cardy 
has shown that the addition of weak disorder on a system of N coupled Ising models ( a 
model presenting a first order transition driven by fluctuations) makes the system flow to 
N decoupled Ising models J/J. This study has then been extended to the case of iV 3— state 
Potts models by Pujol; the result was a non-Ising like second order transition (for N > 2) 
||. The g— states Potts model with q > 4 appears as a more natural model to test and 
analyse more deeply the effects of disorder, because the first order transition is now of mean 
field type, so much stronger. Nevertheless, there is no analytical approach able to control 
the effect of disorder on such systems. On the other hand, Chen et al. have investigated 
the 8-states potts Model using Monte Carlo simulations and confirmed the transition to be 
continuous but also found numerical values of the critical exponents consistent with those 
of the pure Ising model |J. Yet, more recent numerical studies of Cardy and Jacobsen |i"0 | 
and of Picco |TT| are in a clear disagreement with the latter conclusion. They found a mag- 
netic exponent f3/v which varies continuously with q. It has also been shown, using finite 
size scaling combined with conformal invariance, that the values of the central charges are 
related to one another by a factorization law c(q = q\ x q 2 ) = c(qi) + c(q 2 ). Therefore, this 
measure is not able to distinguish between a non trivial random behavior for a q = q\ x q 2 
state random Potts model and two decoupled q\ and q 2 state Potts models. This could sug- 
gest some links between a disordered q = q\ x g 2 -state Potts model and the corresponding 
disordered gi,g 2 -state potts models despite the magnetic exponents appear different for the 
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Up to now, the only analytical results concern several g-state Potts models with disorder. We 
intend in this letter to analyse the case of two different Potts models with disorder. There- 
fore, we investigate analytically (in a perturbative sheme in powers of (q — 2)) the behavior 
of two coupled qi and g 2 -state potts models with q ly q 2 G {2, 3, 4} in order to compare it to 
the numerical results of the q = q\ x g 2 -state Potts model. We find a rather complex situa- 
tion: When two different models are considered, there is a factorisation law and the models 
decouple. Therefore, our analysis confirms no apparent relations (at least perturbatively!) 
between the q = qi x g 2 -state Potts model and the corresponding disordered q±,q 2 state 
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Potts models, despite their central charge to be similar. Nevertheless, when q\ = q2 = 3 or 
4, the system flows in a strong coupling regime, despite the presence of disorder, what is in 
apparent contradiction with the Imry-Wortis argument. 



Our model consists of one gi-state Potts model coupled to a q 2 — state Potts model by 
their energy operators (g, e {2, 3, 4}). The Hamiltonian of the system has the following form 

Hi, H2 are respectively the Hamiltonians of the pure qi and qi state Potts models, m; the 
reduced temperatures (i = 1,2), corresponds to the energy operators of the pure models 
and finally gi2 is the coupling constant associated to the interaction term. The partition 
function can hence be written as 

£ = jy e [Si+S°+gi2 J d?x £ 1 e 2 +m 1 J d?x e 1 +m 2 J d?x e 2 ] /^\ 

where Si^ is a spin operator of the gi-state Potts model. Therefore, a correlation function 
< O(0) O(R) >, where O is some local operator, is expanded perturbatively like: 

< O(0)O(R) >=< O(0)O(R) > + < SiO(0)O(R) > +- < S]O(0)O(R) > + ■ ■ • 

where <>o means the expectation value taken with respect to S® + S® and Si = J d 2 x £162- 
The calculation of correlation functions can be performed with the Coulomb-gas representa- 
tion |12|]. The central charge of a g^-state Potts model is written as c = | + e,, where can 
also be used as a short distance regulator for the integrals involved in correlation functions 
calculations. Then the limit e — > corresponds to the Ising model while the Potts model is 
obtained for some finite value of e. The recent numerical results of Cardy et al. [|K| have 
proved that the (q — 2) expansion gives accurate results for q = 3 and good qualitative results 
for q = 4 (with an error around 4% on the magnetic exponent of the disordered 4— state 
potts model). 

When we have only one coupling constant go, its renormalisation is determined directly by 
a perturbative computation, g is also given by the Operator Algebra (O.A.) producing 
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go 



e a {z)e b {z)d 2 z + 7; (#0 J e a {z)e b (z)d 2 z^j H = g J e a (z)e b {z)d 2 z , (3) 



with g = go + ^2^0 + ' ' ' where A2 comes from the contraction 

- J e a {z)e b (z)d 2 z J e c {z)e d {z)d 2 z = A 2 J e a {z)e b {z)d 2 z (4) 

Therefore, in the limit — » 0, the 2-loop renormalisation group equation associated to (|I|) 
reads 
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Here, q = 2 — A £i , with A £i the physical dimension of the energy operator £j. For more 
details, we refer to references |T2] , ||. In our notation, we have e = for the Ising model, 
e = | for the 3— state Potts model and e = 1 for the 4— state Potts model. From the equation 
we clearly see that the system is driven in a strong coupling regime indicating probably 
a first order transition. A similar situation has already been encountered, when one couples 
two g-state Potts models |§. Moreover, it can be proved exactly that we have a mass gap 
generation jn|. The equation ([5]) can be generalized for one Ising model coupled to several 



Potts models giving a new interesting fixed point structure ]14 



We now add quenched randomness coupled to the local energy densities. This can be done 
by introducing in (JI|) position dependent random mass terms — > mi(x) with rrii(x) = 
and mi(x)mj(y) = Aij5(x — y). Ajj represents the 2x2 symmetric covariance matrix whose 
elements are strictly positive as it should be. If we have considered a diagonal covariance 
matrix corresponding to independent disorders for each models, then A12 would have been 
generated by the R.G. equations. We then apply the replicated method by introducing n 
copies of the system and averaging Gaussian distributions for rrii(x). We finally obtain the 
sum of n q\ and g 2 -state Potts models coupled by their energy densities. The Hamiltonian 
can be written as 



H = H! + H 2 - g 12 J d 2 x Y, £ i £ 2~ J d 2 * 

a <ij,a,()> 



A--e a e b 



(6) 



where i, j = 1, 2 and a, b runs from 1 to n and < . . . > means that when i = j then a 7^ b. 
The hamiltonian (|6|) has four coupling constants. It is a generalization of the study of two 
Ising models with disorder @ or two Potts models with disorder J8|. We now derive the 
2— loops R.G. equations associated to (^). The procedure employed in a generalisation of 
The methods we use are explained and detailed in ref. @]. Nevertheless, since we have 
to compute correlations functions involving different energy operators, new integrals appear 
at one and two loops. Details of the algebra will be presented elsewhere Therefore, 
taking the replica limit n — > 0, the four beta functions are, at two loops 
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d\ogl 

dAu_ 

dlogl 

dA 22 

dlogl 

dA 12 

dlogl 



€j ^ 1 9i2 - 0i 2 (Aii + A 22 ) + -g X2 {Al + A\) + 3A^ 12 + A 12 ^ 2 2 

eiAu - 2A 2 n + 2g l2 A l2 + 2A^ - 2A n A 12 # 12 - A n g 2 12 

e 2 A 22 - 2A 2 22 + 2# 12 A 12 + 2A^ 2 - 2A 22 A 12 # 12 - A 22 ^? 2 (7) 
ei + £2 
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: A 12 + g 12 (A n + A 22 ) - A 12 (A n + A 22 ] 



-(A 2 U + Aj 2 ) + A^ 2 - 2A 12 £ 12 - g\ 2 
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In these equations we have made the change gu — * ^gvi and A — > AttA. We have not 
considered replica symmetry breaking. Moreover, recent analytical and numerical results on 
disordered Potts models are in favor of a non replica symmetry breaking scenario [TO. When 
we couple two similar models (Ising or Potts), we have t\ = e 2 and gu = g, An = A 12 = 
A 2 2 = A and the flow reduces to 

d 9 n.A , J A , ^AJ 



dlogl 



eg - 2gA + g z A + 4gA' 



' lA eA - 2A 2 + 2A# + 2A 3 - g 2 A - 2A 2 g 



dlogl 



We recover in this case the results of Pujol ||. 

We are looking for the fixed points by taking = = 0. For the system there 
are three fixed points 

9 = ; A = (9) 



the trivial one, and 



9 = ; A = ^ + ^ + 0(e 3 ) (10) 
the fixed point corresponding to a decoupling of the models and 

2 2 

9 = j + 0(e 3 ) ; A = i±^ + 0(e 3 ) (11) 
a new non-trivial one mixing both models. 

When we consider one Ising model coupled to a 3 or 4-state Potts model, we take the limit 
ei — > and e 2 = e, we find the following fixed points 

g* 12 = ; AJ 1 =A^ = AJ 2 = (12) 
g\ 2 = 0; A* n = 0; A* 12 = ~ + L + (e 2 ); A^ 2 = + o(e) (13) 

& = 0; A n = 0; A^ 2 = 1 + ^ + o(e 2 ); AJ 2 = 1 + o{e) (14) 



We have also three fixed points, the trivial one flT2|), the one corresponding to a perfect de- 
coupling of the disordered Ising and Potts models ( |13"D and a new fixed point mixing a priori 
both models. Let us notice that A* 2 is undetermined at one loop and two loop calculations 
just enable to compute the first order in e. And finally, when one considers a 3— state coupled 



to a 4— state Potts model, we only find the fixed points ([I4) , fll3Q . 

In order to study the stability of each of these fixed points, we re-express the systems 
(0) around the above solutions. Therefore, we write gi 2 = g* 2 + Sg^, Ay = A*j + SAij 
and keep only the smallest order in e. We thus obtain a linear system 5X = AX with 
X = (5 An, <5A 22 , <5A 12 , 5<7i 2 ). All the information concerning the stability is contained in 
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the matrix A. For the R.G. equations (pj), A is only a 2 x 2 matrix. By calculating the 



eigenvalues of A for each fixed points, we can show that the fixed points (Q), (10) are unsta- 
ble. For the fixed point (|TT|), a numerical study of the flow (fourth order are needed for an 
analytical study of the stability of this fixed point!) shows that it is also unstable contrary to 
what was stated in ||. Consequently, the system is driven in a strong coupling regime. Note 
that this result is surprising because it is in contradiction with Imry and Wortis arguments 
0. Nevertheless, we must not forget that we have done a perturbative analysis at two loops, 
so another non-perturbative fixed point can not be ruled out (or a perturbative one which 
needs three loops terms to be stabilized). Therefore, it would be very interesting to study 
the critical behavior of this model numerically. 

A similar analysis can be repeated for ei = 0. Thus, we can see that (0) is unstable, that 
(|T3|) is stable, and finally that ( |H]) is stable only in the space = 0. In fact, the fixed point 
(|14]) appears as a generalisation of (|T0"|) for different minimal models. We have represented in 
Figure 1, the projection on the flow in the (g\2, A 12 ) plane for two different initial conditions 
(points A and B ). We clearly see that the flow first try to go away and then is driven 
by disorder at the origin corresponding to a perfect decoupling of the models. As it has 
been already noticed in J7|, such a flow is unusual because it violates the c— theorem JL7|. 
There are no corrections to the critical exponents which remain those of the corresponding 
decoupled disordered models. Note that in || , a new stable fixed point was found (for N > 2 
3— state Potts models), but with similar critical exponents that one disordered Potts model. 

In this paper, we have analysed the behavior of two different coupled minimal models 
in presence of disorder. This constitutes a new example of first order transition driven by 
randomness in a second order one (see Figure 1). Moreover, we have shown than the models 



factorize, namely decouple. If we compare the critical exponent (3/v ~ 0.142 found in ]10 
for the disordered 6— state potts model, it appears clearly different from the one of the Ising 
or disordered 3— state Potts model. Therefore, there are no relations at a perturbative level 
between both models despite their central charges are similar. On the other hand, it would 
be very interesting to test numerically the behavior of two coupled qi, g 2 Potts models with 
disorder and to analyse the non-perturbative area in relation with the q\ x g 2 -state Potts 
model. Indeed, for large g\2 , we could imagine a cross-over phenomena in the system of two 
coupled minimal models. The case of two similar Potts models is special because it seems 
to violate the Imry- Wortis argument perturbatively. It deserves a more accurate analysis, 
either analytically from, for example, the integrable point of view or numerically. Finally, 
such an analysis can be extended to the case of N Ising models to M Potts models. We 
expect similar behaviors except, maybe, for the cases N or M = 2. 
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FIGURE CAPTION 



FIG.l : The projection of two flows in the (gn, A i2 ) plane for one Ising model coupled to 
a 3— state Potts models with disorder. A$, Bq correspond to two different initial conditions. 
We clearly see that the flow is driven by randomness at the origin. It corresponds to a 
decoupling regime. 
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